Nonlinear Modes of a Macroscopic Quantum Oscillator 
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We consider the Bose-Einstein condensate in a parabolic trap as a macroscopic quantum oscillator 
and describe, analytically and numerically, its collective modes-a nonlinear generalisation of the 
(symmetric and antisymmetric) Hermite-Gauss eigenmodes of a harmonic quantum oscillator. 
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The recent observation of different types of Bose- 
Einstein condensation (BEC) in atomic clouds Q led to 
the foundation of a new direction in the study of macro- 
scopic quantum phenomena. From a general point of 
view, the dynamics of gases of cooled atoms confined in a 
magnetic trap at very low temperatures can be described 
by an effective equation for the condensate wave function 
known as the Gross-Pitaevskii (GP) equation This 
is a classical nonlinear equation that takes into account 
the effects of the particle interaction through an effective 
mean field, and therefore it can be treated as a nonlinear 
generalization of a text-book problem of quantum me- 
chanics, i.e. as a macroscopic quantum oscillator. 

Similar models of the confined dynamics of macro- 
scopic quantum systems appear in other fields, e.g. in 
the case of an electron gas confined in a quantum well 
||, or optical modes in a photonic microcavity Q. In 
all such systems, confined single-particle states are re- 
stricted to a set of discrete energies that form a set of 
eigenmodes. A classical and probably most familiar ex- 
ample of such a system is a harmonic quantum oscillator 
with equally spaced energy levels j^] . 

When, instead of single-particle states, we describe 
quasiclassically a system of interacting bosons in a 
macroscopic ground state confined in an external poten- 
tial, a standard application of the mean-field theory al- 
lows us to introduce a macroscopic wave function as a 
classical field ^(R, t) having the meaning of the order 
parameter. The equation for the function ^(R, t) looks 
similar to that of a single-particle oscillator, but it also 
includes the effect of interparticle interaction, taken into 
account as a mean-field nonlinear term. Then, the im- 
portant questions are: Does the physical picture of eigen- 
modes remain valid in the nonlinear case, and what is the 
effect of nonlinearity on the modes ? In this paper we 
analyse nonlinear eigenmodes of a macroscopic quantum 
oscillator as a set of nonlinear stationary states that ex- 
tend the well-known Hcrmite-Gauss eigenfunctions. Wc 
also make a link between seemingly different approxima- 
tions, the well-known Thomas- Fermi approximation and 
the perturbation theory developed here for the case of 
weak nonlinearity. For both attractive and repulsive in- 
teraction, we demonstrate a close connection between the 



nonlinear modes and (bright and dark) multi-soliton sta- 
tionary states. 

We consider the macroscopic dynamics of con- 
densed atomic clouds in a three-dimensional, strongly 
anisotropic, external parabolic potential created by a 
magnetic trap. The BEC collective dynamics can be de- 
scribed by the GP equation, 
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where ^(R, t) is the macroscopic wave function of a con- 
densate, V(R) is a parabolic trapping potential, and 
the parameter Uq — AirTi 2 (a/rn) characterises the two- 
particle interaction proportional to the s-wave scattering 
length a. When a > 0, the interaction between the parti- 
cles in the condensate is repulsive, whereas for a < the 
interaction is attractive. In fact, the scattering length a 
can be continuously detuned from positive to negative 
values by varying the external magnetic field near the 
so-called Feshbach resonances || . 

First of all, we derive from Eq. ([l]) an effective 
one-dimensional model, assuming the case of a highly 
anisotropic (cigar-shaped) trap of the axial symmetry 
V(R) = \moj\ {R\ + XX 2 ), where R ± = VY 2 +Z 2 . 
This means that A = wjj/o;^ <C I, and the transverse 
structure of the condensate, being close to a Gaussian in 
shape, is mostly defined by the trapping potential 0. 

Measuring the spatial variables in the units of 
the longitudinal harmonic oscillator length a/ lo = 
(ft/mwvA) 1 / 2 , and the wavefunction amplitude, in units 
of (Sw_l/2[/oVA) 1 ^ 2 , we obtain the following dimension- 
less equation: 
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V 2 * - [X^{y 2 + z 2 ) + x 2 }® + ct|*| 2 ^ = 0, (2) 



where time is measured in the units of (2/uij_y/\), 
(x,y,z) = (X,Y, Z)/aho, and the sign a = sgn(et) = ±1 
in front of the nonlinear term is defined by the sign of 
the s-wave scattering length of the two-body interaction. 

We assume that in Eq. (Q) the nonlinear interaction is 
weak relative to the trapping potential force in the trans- 
verse dimensions, i.e. A« 1. Then, it follows from Eq. 
(0) that the transverse structure of the condensate is of 
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order of A, and the condensate has a cigar-like shape. 
Therefore, we can look for solutions of Eq. (||) in the 
form, 

V(r,x,t) = <P{r)ijj(x,t)e~ 2i ^\ 



where r = y^y 2 + z 2 , and <f>(r) is a solution of the aux- 
iliary problem for the 2D radially symmetric quantum 
harmonic oscillator 

Vi$ + 2 7 $ - (r 2 /A)$ = 0, 

which we take in the form of the no-node ground state, 
$ (r) = Cexp(-7r 2 /2), where 7 = 1/V% To preserve 
all the information about the structure of the 3D conden- 
sate in an asymmetric trap describing its properties by 
the effective GP equation for the longitudinal profile, we 
impose the normalisation for $o(r) that yields C 2 = 'y/n. 

After substituting such a factorized solution into Eq. 
(0), dividing by $ and integrating over the transverse 
cross-section of the cigar-shaped condensate, we finally 
derive the following ID nonstationary GP equation 



if + 0-^ + ^ = 0. 



(3) 



The number of the condensate particles N is now de- 
fined as N = (hui/2UoVX)Q, where 
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\ip(x,t)\ 2 dx 



(4) 



is the integral of motion for the normalised nonstationary 
GP equation (|). 

Equation (|3) includes all the terms of the same order, 
and it describes a longitudinal profile of the condensate 
state in a highly anisotropic trap. In the linear limit, 
i.e. when formally a — > 0, Eq. (||) becomes the well- 
known equation for a harmonic quantum oscillator. Its 
stationary localised solutions, 



ip(x, t) — <f>(x)e 
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exist only for discrete values of fl, such that 

0„ = l + 2n, n= 0,1,2,..., 

and they are defined through the Hermite-Gauss 
polynomials, (f) n (x) — c n e~ x ' 2 H n (x), where c„ = 
(2"n! y^)- 1/2 , and 



2 , d n (?- x2 / 2 ) 
H n (x) = (-l)" c * 2 /2li£ >_ 

K ' K ' dx n 



(6) 



so that H = 1, Hi = 2x, etc. 

In general, the localised solutions of Eq. (|^) for a 7^ 
can be found only numerically. All such solutions can be 
characterised by the dependence of the invariant (EJ) on 
the effective nonlinear frequency Q. However, in some 



particular limits we can employ different approximate 
methods to find the localised solutions analytically. 

First of all, to describe the effect of weak nonlinearity, 
we use the perturbation theory based on the expansion 
of the general solution of Eq. (pi) in the infinite set of the 
eigenf unctions (||). A similar approach has been used 
earlier in the theory of the dispersion-managed optical 
solitons ||. To apply such a perturbation theory, we 
look for solutions of Eq. (||) in the form [cf. Eq. (ph] 



ip(x,t) 
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where <p n (x) are the eigenfunctions of the linear equation 
for a harmonic oscillator that satisfy the equation 



dx 2 



x 2 (j) n + Q n (/>n = 0. 



(8) 



Inserting the expansion (0) into Eq. (||), multiplying by 
(j) n and averaging, we obtain a system of algebraic equa- 
tions for the coefficients 



(CI — fl m )B m — a V m . n j t kB n BiBk — 0, 
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where 
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4>m {x)4> n (x)(pi {x)(j) k (x)dx. 



Equation (^) can be also rewritten in the traditional form 
S(H + QQ) = 0, and it allows us to develop a pertur- 
bation theory for small nonlinearities, for any Hermite- 
Gaussian eigenmode. For example, let us consider the 
ground state mode at n = 0. Assuming Bq B m 
for m ^ and the condition of the symmetric solution 
4>{x) — x), i.e. i?2fe+i = for any k, we find the 
corrections, fl w flo + erVb,o,o,o|-6o| 2 , and 

B 2k = £^\B \ 2 B , k^0. 

This allows us to calculate the asymptotic expansion of 
the invariant Q for small nonlinearities, 

Q = J2 l B2fe l 2 ~ -™o(« - n )[i + b (n - a,) 2 ], (10) 



where the coefficients are: 
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V2^, b = J2 



(2k)! 



^ (4k) 2 (k\2 2k ) 2 



Higher-order modes can be considered in a similar way, 
and the results are similar to Eq. (^) where do and 60 
change to a n and b n respectively, where a n and b n depend 
on the mode order. 

In the opposite limit, i.e. when the nonlinear term or 
potential are large in comparison with the kinetic term 
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given by the second-order derivative, we can use two 
different approximations for describing localized modes. 
For a = +1 and large negative fl, localised modes are 
described by the stationary solutions of the nonlinear 
Schrodinger (NLS) equation, that appears when we ne- 
glect the trapping potential. The NLS one-soliton solu- 
tion is 4> s {x) = a/ — 20 sech(a;-\/— £1), so that the depen- 
dence Q(fi) coincides with the soliton invariant Q s = 
4V— For a — — 1 and large positive J7, the ground- 
state solution can be obtained by using the so-called 
Thomas-Fermi approximation, based on neglecting the 
kinetic term - this yields 4>tf (x) ~ — x 2 . 

In general, we should solve Eq. (||) numerically. Fig- 
ures 1(a) to 1(b) present examples of the numerically 
found ground-state solutions of Eq. (g) in the form (||) as 
continuous functions of the dimensionless parameter O, 
for both negative (a) and positive (b) scattering length. 
For — ► 1, i.e. in the limit of the harmonic oscillator 
ground-state mode, the solution is close to Gaussian for 
both the cases. When deviates from 1, the solution pro- 
file is defined by the type of nonlinearity. For attraction 
(a = +1), the profile approaches the sech-type soliton, 
whereas for repulsion (a = —1) the solution flattens, and 
it is better described by the Thomas-Fermi approxima- 
tion, that is good except at the edge points. 




FIG. 1. Examples of the condensate density \ip\ 2 for the 
first three nonlinear modes described by the stationary solu- 
tions of Eq. (^) for the negative (a = +1, upper row) and 
positive (a = —1, lower row) scattering length, respectively. 
The values of Q are given next to the curves. 

In Fig. 2 we present the dependence of the invariant 
Q on the parameter f2, for both the types of the ground- 
state solution, corresponding to two different signs of the 
scattering length. The dashed curve for the zero-order 
mode (marked as Oth-mode) shows the dependence Q S (Q) 
for the soliton solution of the NLS equation without a 
trapping potential. In the asymptotic region of negative 



f2, i.e. say for il < —2, the dependence Q(fi) for the 
BEC condensate in a trap approaches the curve Q s (£l). 
This means that for such a narrow localised state the 
effect of a parabolic potential is negligible, and the con- 
densate ground state becomes localised mostly due to an 
attractive interparticle interaction. In contrast, for large 
positive the effect of a trapping potential is crucial, 
and the solution of the Thomas-Fermi approximation, 
<^TF(aO = — x 2 , defines the common asymptotics for 
all the modes, Qtf ~ |i7 3 ^ 2 (dashed-dotted curves in 
Fig. 2). 
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FIG. 2. Invariant Q vs. Q for the first three nonlinear 
modes in the case of attraction (a = +1, left), and repul- 
sion (a = —1, right). Dashed curves for a = +1 are de- 
fined by the soliton invariants Q s , 2Q S , and 3Q S , respectively. 
Dashed-dotted curves for a = — 1 are the asymptotic limit 
given by the Thomas- Fermi approximation, Qtf = |f2 3,/2 . 
Dotted lines show the results of the perturbation theory ap- 
plied to the ground state, see Eq. (px|). 

As has been mentioned above, in the linear limit 
(a —>■ 0), Eq. (||) possesses a discrete set of localised 
modes described by the Hermite-Gauss polynomials. We 
have demonstrated that all such modes can be readily 
calculated by the perturbation theory in the weakly non- 
linear approximation, and therefore they should all exist 
for the nonlinear problem as well, describing an analyti- 
cal continuation of the Hermite- Gauss linear modes to a 
set of nonlinear stationary states. In application to the 
BEC theory, these non-ground-state solutions were first 
discussed by Yukalov et al Q. 

Figures 1(c) to 1(f) show examples of the first- and 
second-order modes for both negative and positive scat- 
tering length, respectively. In the limit Vl — > fl n , all those 
modes transform into the corresponding eigenfunctions of 
a linear harmonic oscillator. 

It is clear that nonlinearity has a different effect for 
the negative and positive scattering length. For the 
negative scattering length (attraction), the higher-order 
modes transform into multi- soliton states consisting of a 



3 



sequence of solitary waves with alternating phases [see 
Figs. 1(c) and 1(e)]. This is further confirmed by the 
analysis of the invariant Q vs. fl, where all the branches 
of the higher-order modes approach asymptotically the 
soliton dependencies Q n ~ (n + l)Q s , where n is the 
order of the mode (n = 0, 1, . . .). From the physical 
point of view, in the case of attractive interaction the 
higher-order stationary modes exist due to a balance be- 
tween repulsion of out-of-phase bright NLS solitons and 
attraction imposed by the trapping potential. The analy- 
sis of the global stability of such higher-order multihump 
multi-soliton modes is still an open problem, however the 
recent results indicate that, at least in some nonlinear 
models, multihump soliton states can be stable JhJ. 

For the positive scattering length (a = —1), the higher- 
order modes transform into a sequence of dark solitons 
(or kinks) jll| , so that the first-order mode corresponds to 
a single dark soliton, the second-order mode, to a pair of 
dark solitons, etc. [see Figs. 1(d) and 1(f)]. Again, these 
stationary solutions satisfy a force balance condition - re- 
pulsion between dark solitons is exactly compensated by 
an attractive force of the trapping potential. 

The modal structure of the condensate macroscopic 
states described above and summarised in Fig. 2 for both 
positive and negative values of the scattering length al- 
lows us to draw an analogy between BEC in a trap and 
the guided-wave optics where the condensate dynamics 
in time corresponds to the stationary mode propagation 
along an optical waveguide, with the parameter f2 as 
the propagation constant. As is well known from dif- 
ferent problems of guided-wave optics, in the presence 
of interaction the guided modes become coupled and the 
coupling can lead to both power exchange and nonlinear 
phase shifting between the modes. In application to the 
BEC theory, the mode coupling resembles a kind of inter- 
nal Josephson effect. These issues are beyond the scope 
of this paper and will be analysed elsewhere Jl2] . 

The theory of nonlinear stationary modes of a macro- 
scopic quantum oscillator developed above for the ID 
analog of BEC can be easily extended to both 2D and 
3D cases. Moreover, the coupled-mode theory for a sin- 
gle condensate is closely connected to the dynamics of 
strongly coupled two-component BECs [fi"3| where exci- 
tation of an antisymmetric (or, in our notation, the first 
order) collective mode, in the form of collapses and re- 
vivals, has been recently observed experimentally 

At last, we would like to mention that the basic con- 



cepts and results presented above can find their applica- 
tions in other fields. For example, the effect of nonlinear- 
ity can lead to a mode coupling in the so-called "photonic 
atom" , a micrometer-sized piece of semiconductor that 
traps photons inside 0, or two such 'photonic atoms' 
coupled together, a 'photonic molecule' |Q. In such pho- 
tonic microcavity structures, the macroscopic nature of 
the states may lead to different nonlinear effects includ- 
ing the mode mixing and power exchange. 

In conclusion, we have analysed nonlinear stationary 
modes of a macroscopic quantum oscillator considering, 
as an example, the cigar-shaped Bose-Einstein conden- 
sate in a parabolic trap. 
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